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Abstract Collisions between particles suspended in a fluid play an important role 
in many physical processes. As an example, collisions of microscopic water droplets in 
clouds are a necessary step in the production of macroscopic raindrops. Collisions of dust 
grains are also conjectured to be important for planet formation in the gas surrounding 
young stars, and also to play a role in the dynamics of sand storms. In these processes, 
collisions are favoured by fast turbulent motions. Here we review recent advances in 
the understanding of collisional aggregation due to turbulence. We discuss the role of 
fractal clustering of particles, and caustic singularities of their velocities. We also discuss 
limitations of the Smoluchowski equation for modelling these processes. These advances 
lead to a semi-quantitative understanding on the influence of turbulence on collision rates, 
and point to deficiencies in the current understanding of rainfall and planet formation. 
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1 Introduction 


Fluids which are encountered in nature often carry a suspension of small particles [T] , which 
could be dust grains, small liquid droplets or objects of more complicated shape, such as 
ice crystals w, or even living organisms |3j . Their presence is typically evident from the 
optical behaviour of the fluid due to scattering of light by the suspended particles, and 
this is the reason why the atmosphere or the sea are often opaque or cloudy. In clouds, the 
absorption of short-wavelength radiation from the sun, and of long-wavelength radiation 
from the earth, play an essential, yet incompletely understood, role in the energy budget 
of the planet [5]. By scavenging of aerosol particles and of other gases, water droplets 
contribute in an essential way to the dynamics of the atmosphere and of the climate. 
Reliably modelling the distribution of particle sizes in suspensions is therefore essential to 
understand these processes. 

As well as influencing optical properties, the stability of suspended particles in large 
bodies of fluid can have important implications. The formation of rain drops by the 
coalescence of vast numbers of microscopic water droplets which make up atmospheric 
clouds, is one of the central questions facing cloud microphysics, which continues to receive 
much attention Also, the standard model for formation of planets depends upon 

the collision and coalescence of dust grains in the atmosphere around a young star [7]. 
Explaining the existence of our planet, and the weather phenomena that make it habitable, 
depends upon the instability of aerosol suspensions to collisions. 

Attempts to make a quantitative theory for rain initiation or planet formation, however, 
run into difficulties. The collision rates of the atmospheric aerosol particles, resulting from 
differential settling rates of water droplets of different sizes or from Brownian motion, 
appear to be insufficient to explain the rapid onset of rain from many types of cloud. 
Similarly, it is difficult to argue that the collision rate of dust grains is adequate to explain 
planet formation. 

One possible route to resolving these problems involves turbulence Turbulent 

motion is a robust phenomenon, which occurs in many situations where a large body of 
fluid is in motion. It is well known that the dispersion of small particles in a turbulent 
environment is much faster than can be achieved by molecular diffusion |10] . which suggests 
that turbulence could also greatly enhance collision rates. This argument implies that 
turbulence could also play an essential role in the coagulation process of particles. 

The mechanisms whereby turbulence can enhance the collision rate have only become 
clear in the last few years. This review will explain the current understanding, which is 
forming a coherent and internally consistent picture, supported by numerical experiments 
on accurate simulations of turbulent flow. The theoretical picture of the collision rate 
proved to be multifaceted, involving concepts from dynamical systems theory, fractal ge¬ 
ometry, stochastic processes and optics, as well as results from fluid dynamics and the 
theory of turbulence. 

While some aspects of quantifying collision rates of small particles in turbulent flows 
may require further work, it seems that the underlying physical principles are now qualita¬ 
tively well understood. The knowledge gained is expected to lead to a deeper understand¬ 
ing of collective behaviours in turbulent suspensions. It remains to apply this knowledge 
to significant problems such as explaining rainfall, planet formation, and properties of 
particle-laden turbulent flows such as sandstorms or powder-snow avalanches [H]- At this 
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Figure 1: Particles in a turbulent flow can show pronounced clustering, which samples a 
fractal measure. The image is a two-dimensional simulation of particles in an incompress¬ 
ible flow. The simulation includes particles with three different masses, shown in different 
colors. 


level, fundamental problems remain outstanding. The general framework originally devel¬ 
oped by Smoluchowski [12] to describe coagulation processes in suspensions of Brownian 
particles, based on a mean-field approximation, may appear as an enticing starting point. 
The rapid increase of the collision rate when the size of the particle increases, as it hap¬ 
pens in the case of settling droplets in a cloud, can lead to runaway growth, which is a 
feature of the formation of raindrops from microscopic water droplets. This phenomenon 
is known as gelation in the polymer physics literature. Surprisingly, it has been shown 
that when gelation is modelled using the Smoluchowski equation, the time required for the 
gelation transition may be strictly equal to zero m- This instantaneous gelation is clearly 
unphysical; it implies that mean-field descriptions based on the Smoluchowski approach 
have to be applied with caution. 

We provide here a critical discussion of results concerning the rate of rain drop and 
planet formation, which shows that while turbulence can dramatically increase the collision 
rate, it is not entirely clear whether it is sufficient to explain the coagulation rates observed 
in nature, and alternative physical explanations may be required. 

Over the past few years, the subject of collisions induced by turbulence has received 
a surge of interest, and as such, it has been the subject of several other review articles, 
focusing mostly on meteorological applications |15l [T6l [6]. This review is more focused 
on the lessons learned from a general fundamental physics perspective. It is our belief 
that some of the results shown here will be relevant in other fields of physics, and thus of 
lasting interest, beyond the specific motivations of the original study. 

This article is organised in the following way. In section [2] we review elementary 
material concerning the motion of particles in a fluid flow, and the collision rate in a 
suspension, independently of the precise nature of the flow. The more specific aspects of 
turbulent flows, necessary for the purpose of this review, are discussed in section [Sj Section 
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Figure 2: a The distribution of particles in a turbulent flow can show singularities 

which result from the projection of folded manifolds in phase space: this is a simulation 
of a two-dimensional compressible flow, taken from HU. b The singularities are termed 
caustics, because they have the same structure as singularities in optics. The image (a 
photograph by Amanda W. Peet) shows a pattern of caustics on the bottom of a swimming 
pool, resulting from partial focussing of sunlight due to curvature of the water surface by 
irregular waves. The two images are remarkably similar. 

lUwill introduce various mechanisms whereby turbulence can influence the collision rate, 
including particle clustering and the effects of caustics forming in the phase space of the 
suspended particles. These diverse contributions to the collision rate are synthesised into 
a unified approximation scheme in section [5l which is shown to provide a very accurate 
description of the results of numerical simulation. Having discussed the determination for 
collision rates we turn to discussion of applications. Section [6] considers the question of 
whether aerosols undergo a gelation transition, and whether the Smoluchowski equation 
provides an appropriate description. Sections [7] and [8] consider applications to rainfall and 
planet formation. In both cases we argue that the insights from considering turbulent 
enhancement of collision rates do not appear to be sufficient to resolve all of the problems 
with understanding these processes. The prospects for further developments are considered 
in section [9l which is our conclusion. A detailed analysis of both fractal clustering and 
caustics, which are crucial to the determination of the collision rate, is provided in the 
case of a solvable one-dimensional model in Appendix A (Section I lOj) . 

2 Definitions and equations 

2.1 Particle motion in a fluid flow 

The determination of the motion of particles, even in simple (laminar) flow configurations, 
is a difficult task. In the applications we have in mind, particles have small sizes, and this 
leads to a solvable problem. In this limit, the flow can be approximated as constant over 
a domain much larger than the particle, and the Reynolds number of the particle is so 
small that the nonlinear term in the Navier-Stokes equations drops out, so the problem 
can be explicitly solved HU IS]- The equations of motion involve several terms, which 
result from the viscous drag of the particle, the gravitational settling, pressure effect, the 
so-called added mass, and an history (Basset-Boussinesq) force. 

Rain drops in a cloud or particles in the interstellar medium have densities, pp, which 
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are much larger than the fluid density, pf. pp/pf S> 1. The resulting inertia of the particles 
may be large enough to prevent them from exactly following the flow. In this limit, it 
has been demonstrated m that the two dominant forces on small particles are due to 
viscous drag, which causes the particle velocity to relax towards that of the fluid, and to 
gravitational settling. The equations of motion for small spherical particles of radius a are 
determined by Stokes’ formula, so that the equation of motion is 


r = V , 


V = — [u{r, t) — v] + g 


( 1 ) 


where 

— - — 82. 

9 v Pi 

is the particle relaxation time, determined from Stokes’ formula for the drag on a moving 
sphere (in ([2]), v is the kinematic viscosity). The equations of motion (|1I2I) are valid only 
in the limit where the suspended particles are very small and very dense: Pp/pi 1. 
When the ratio Pp/pi is smaller than ~ 10, the history force becomes important [20]. In 
cases where a droplet moves through a fluid which has a comparable viscosity, equation 
([^ must be modihed 1211122]. 

In astrophysical applications, the gas phase may have extremely low density, so that 
the size of the aerosol particles may be much smaller than the mean free path of the 
gas. In this case equation ([T]) still provides an accurate description of the motion, but the 
relaxation time Tp is given by a different expression: 


Tp = Ke 


CsPg 


( 3 ) 


where Cg is the velocity of sound in the gas with density pg, and Ke is a dimensionless 
constant which is of order unity [23] . The value of Ke depends on assumptions about the 
energy transferred when a gas molecule collides with the surface of the particle. 

Describing the motion of larger particles in a turbulent flow remains an outstanding 
challenge, whose numerical and experimental study is at an early stage [2ll[25l[26l[271[28]. 


2.2 Collision rates 

In principle, the calculation of the collision rate appears to be a straightforward exercise 
in elementary kinetic theory [29], but in practice it is a surprisingly complex problem. 
Before discussing the details, it is necessary to consider some fundamental definitions. 

Our objective is to quantify the rate R for collision of a given particle with any other 
particle in the suspension. This quantity is expected to be a function of the radius a of the 
particle concerned, so our objective is to calculate R{a), which has dimensions of inverse 
time: [R] = T“^. The rate of collision is expected to be proportional to the number 
density of particles in the suspension. Let N{a) 6a be the number density of spherical 
particles with radius in the small interval [a, a + da]. The rate of collisions may be written 

poo 

R{a)= / da'r(a, o') iV(a') (4) 

Jo 

where r(a, a') is a characteristic property of the fluid motion, which is termed the colli¬ 
sion kernel. Thus, the specific problem of determining the collision rate R{a) in a given 
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suspension is solved by determining the collision kernel r(a, a'), which is independent of 
the particle density, and then applying equation (jj]). 

The collision kernel may be expressed in terms of other variables, for example when 
we formulate the Smoluchowski equations in section [6] it will be useful to represent the 
particle size distribution in term of masses, so that the number density of particles in the 
interval [m, m + 5m] is N{m) 5m. In this case equation Q is replaced by 

roo 

R{m) = / dm' K{m,m')N{m') (5) 

Jo 

where K{m,m') is collision kernel expressed in terms of mass. 

In a monodisperse suspension, where the particles all have approximately the same 
radius a, and where their number density is n, equation Q reduces to: 

R = r(a, a) n (6) 

where 

POO POO 

n= da N{a) = / dm N{m) (7) 

Jo Jo 

is the number density of particles. 

The collision rate between particles, which are assumed to be spherical, is the rate at 
which the separation vector between the centres of the particles cross a sphere of radius 
01 + 02 , as illustrated in figured Thus, the collision kernel is expected to be proportional to 
the area of the spherical surface, 47r(oi + 02 )^. It is also proportional to a suitably defined 
average of the relative velocity of particles on the sphere of radius oi + 02 , which will be 
denoted by (|Au|). Also, if the particles have a tendency to cluster together, the collision 
kernel would be proportional to the two-point correlation function C{r) for particles with 
separation r = oi + 02 . Putting these factors together, the collision kernel can be written 
in the form 

Pi, 2 = 2 + a2fC{ai + 02 ) (|Aul) . (8) 

The factor ^ in Q is included because only half of the particles crossing the surface are 
travelling inwards, but this is really just a matter of convention in the definition of (|Aul). 
With a suitable definition of (|Au|), this can be presented as an exact equation, see [301I3T] . 

The collision kernel in (l8|) depends on the mean relative velocity of the particles, (| Au|), 
which is difficult to calculate. The complex problem of determining the relative velocity 
between small particles is simply tractable in two limits, relevant to the problem studied 
here. One case is where the relative velocities might be approximately independent of 
the separation of the particles, as illustrated in figure |3|^a) . This is only possible when 
the particles are able to move relative to the fluid. The other limiting case is where the 
particles are advected with the fluid. For very small advected particles, their relative 
velocity at collision Av is proportional to their separation Ax: 

Av = A Ax (9) 

where A is the matrix of velocity gradients, with elements Aij = 

The approach discussed here is purely geometric, and misses important effects. When 
two particles get close to each other, the fluid trapped between them leads to a lubrication 
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Figure 3: The rate of collision between particles of radius oi and 02 is determined by 
integrating the relative velocity over a spherical surface of radius ai + 02 , indicated by a 
dotted line. 



Figure 4: The relative velocity of particles might be independent of their separation, 
illnstrated in (a) for a larger particle overtaking a smaller one as they both fall under 
gravity. Alternatively, the relative velocity may be proportional to their velocity gradient; 
(b) shows a collision induced by a simple shear flow (in both cases the arrows indicate the 
velocity of the smaller particle relative to the larger one). 
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film, which potentially significantly reduce the collision rates. This reduction can be taken 
into account by introducing a collision efficiency [21 [32]. It will be discussed in Section 16.11 
In the following, we discuss three useful examples of physical situations, amenable to 
an explicit determination of the collision kernel. 

2.2.1 Collision in a gas of particles with a Gaussian distribution of velocity 

For future reference, we estimate here the collision rate between particles spatially uni¬ 
formly distributed in a fluid, with statistically independent Gaussian (Maxwellian) distri¬ 
bution of velocities, such that {v) = 0 and a variance {v'^)a, which depends a-priori on 
the radius a of the particles. The assumption that particles are uniformly distributed 
ensures that C(ai - 1 - 02 ) = 1. Assuming that the distribution of the velocity for particle i is 
P{vi) = ( 27 r(u^) . 6xp[—3ui^/(2(u^)aJ], an elementary calculation shows that the aver¬ 
age velocity required in ([8|) is (|Aur|) = [f )ji/ 2 _ i 0 ads to the following 

collision rate in a suspensions of particles with a Gaussian distribution of velocities: 

ri,2 = (f )''"(<.! + + {t-")™!'-'" (10) 

consistent with classical estimates [33]. 


2.2.2 Collision between settling droplets in still air 

We consider now the case of polydisperse suspension of water droplets settling in still air. 
Integrating ([I]), and taking into account the dependence on the size a of Tp leads to the 
following equations for the settling velocity 


Vs 


-VsCz, Vs = Ka 


= IPRi. 

9 pgu 


( 11 ) 


This assumes that the Reynolds number based upon the particle size is small, which is 
valid for the microscopic water droplets in clouds [32]. Equation m shows a strong 
dependence of Vg on a: large particles settle much faster than smaller ones. Consider the 
collision rate between particles of sizes oi and 02 . In this case, the value of C(ai -|- 02 ) 
reduces to 1, and a simple calculation leads to (|Aur-|) = «:|of — 021 / 2 , which immediately 
leads to: 

ri ,2 = 7 rK(ai-I-02)^101 - a|| . ( 12 ) 

This expression, which grows as a power oc af, when oi 02 , i.e., with a power 4/3 of 
the volume of the larger particles, implies that the collision rate of large particles grows 
very rapidly as their size increases. This fast growth will have an important consequence 
when studying the problem of coagulation, see Section [Gj 


2.2.3 Collision in a simple shear flow 

Whereas in the two previous examples, the motion of the fluid was not playing any role, 
consider now a suspension of particles simply transported in a simple (laminar) shear 
flow: u = S(y,0,0), so particles are moving in the x-direction, with a velocity v = u, 
which depends only on its y-component (see flgureH]). This problem was flrst treated by 
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Figure 5: The velocity gradient is not usually a simple shear. A locally hyperbolic relative 
velocity also induces collisions. 


Smoluchowski [12]. Once again, assuming a uniform distribution of particles in the flow 
implies that C{ai + 02 ) = 1. An elementary estimate of the Avr leads to the following 
expression for (|Aur|) = ||5|(ai + 02 ), which leads to the collision rate: 

ri,2 = ^|5|(ai+a2f . (13) 

The collision rate, which is proportional to the shear rate, S, and the distance (oi + 02 ) 
to the third power, is characteristic of the collision rates obtained as a result of shearing 
motion, at least when the particles follow the flow. In this case, the simple flow structure 
appears through the shear rate, S. 

In turbulent flows, velocity differences are created at many scales, generating, among 
other things, strong shearing motions on small scales. The velocity gradient in turbulence 
or other complex flow field is described by a real-valued 3x3 matrix A, which is generically 
not a simple shear. Other possibilities include the case of locally hyperbolic flow, as 
illustrated in figure O To gain further insight on the statistics of (|Au|) requires a more 
precise description of the properties of the turbulent velocity field, which will be considered 
in the next section. 


3 Elementary properties of turbulent flows 

The enhancement of the collision rate by a simple flow, such as a shear flow considered in 
the previous subsection, provides the first hint concerning enhancement of collision rates 
by turbulent flows. In fact, significant velocity gradients are ubiquitous in turbulent flows. 
This section is aimed at providing an elementary discussion of the properties of turbulence 
relevant to the present discussion. 

In this review, we are concerned with turbulence in an incompressible, Newtonian fluid. 
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described by the Navier-Stokes equations: 


(14) 

(15) 


dt u + (u ■ V)u =——Vp + + f 

Pf 

V- u = 0 . 

where p is pressure. The viscous term dissipates kinetic energy in the fluid, so in the 
absence of any forcing, / = 0, the motion simply decays. The Reynolds number, deflned 
as Re = ULjv, where U and L are the velocity and length scales at which the fluid is forced, 
measures the ratio between the nonlinear term, {u-\/)u and the viscous (dissipative) term, 
In turbulent flows, the Reynolds number is effectively a measure of the intensity of 
turbulence: turbulence is expected to occur whenever the Reynolds number is very large: 
Re S> 1. Thus, the statement Re S> 1 effectively means that, at the forcing scale, viscous 
dissipation does not play much role. Primarily because the kinematic viscosity is a small 
quantity, very large Reynolds numbers are common: for example, a convective instability 
of the atmosphere can create a flow with Re > 10^. 

3.1 Scales in a turbulent flow 

Turbulence is a notoriously difficult problem |34] but we argue that most of what we need 
to know can be surmised from dimensional considerations. Kolmogorov |351l36j introduced 
the powerful notion that the small structures of the flow have very little ‘memory’ of how 
the flow was generated. He argued that only one quantity, the rate of kinetic energy 
dissipation per unit mass, e, is required to characterise fully developed steady and homo¬ 
geneous turbulent motion over a wide range of length scales (termed the ‘inertial range’ 
in the turbulence literature). The quantity e, the power dissipated per unit mass, and 
designated in the following in short as ‘energy dissipation’, is determined by the manner 
in which the turbulence is generated, by means of large scale flows with characteristic 
length scale L and velocity scale U. Dimensional considerations imply that e U^/L. 
In the inertial range of length scales, where e is the only relevant parameter, statistics of 
the flow can be determined by dimensional analysis. For example, in the inertial range 
of separations, the variance of the velocity difference between two points depends only 
upon the separation R = |i?| and the kinematic viscosity ly. Dimensional consistency then 
implies that 

{[u{R,t)-u{0,t)f)=C{e\R\f^ (16) 

where the brackets in (|16l) refer to an average over many flow realizations, and (7 is a 
universal dimensionless coefficient. 

The turbulence generates successively finer scale eddies until the structures become so 
small that gradients increase, making the power dissipated per unit mass, e ~ nu ■ 
significant. The smallest scale reached by the flow is the ‘Kolmogorov lengthscale’, p, 
with a characteristic time scale known as the Kolmogorov timescale, tk- These quantities 
depend only on e and n. Dimensional considerations then imply that 

p = (i/^/e)^/'^ , tk = . (17) 

As finer scales are generated by the flow, it is expected that the statistical properties 
of the flow become homogeneous and isotropic [351 [36] . In this case the kinetic energy 
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dissipation can be obtained directly from (ITO : 

^ = ( 18 ) 
hj 

where Aij, the velocity gradient, is defined by Q. This implies that the typical size of 
the velocity gradient is the inverse of the Kolmogorov time: dujdx ~ . 

3.2 Velocity gradient statistics 

In turbulent flows, the relative motion of two small particles approaching (colliding with) 
each other is ultimately dominated by the velocity gradient tensor. A, as explained in 
Subsection ESI For this reason, we briefly discuss some elementary statistical properties 
of the velocity gradient. 

The isotropy of the flow imposes that the tensor {Aij{x)Aki{x)) is expressible in terms 
of Kronecker 6 tensors. Using the incompressibility condition tr(A) = 0, as well as relation 
(USD leads to the following expression for {Aij(x)Ai^i{x)): 

{Aij (^x'jAf^li^x')') — ^il^jk 

Such an estimate is crucial in establishing elementary results, such as the Saffman-Turner 
collision rate, Eq. (I25|). 

As a particle is transported by the flow, the strain and the vorticity along its trajectory 
decorrelate with a correlation time of the order of tk [371 [38] . 

3.3 The Stokes number 

Comparing the Kolmogorov time tk with the response time of the particles, Tp, provides 
a way to quantify the effect of inertia. This motivates the definition of the Stokes number: 

St = . (20) 

Tk 

For St <C 1, particles are advected by the fluid, and collisions are the result of the 
shear (the relative motion). When St S> 1, the inertia of the particles allows them to 
move relative to the surrounding fluid, thus leading to entirely different phenomena. The 
Stokes number is the single dimensionless parameter which distinguishes different physical 
regimes of the collision process. 

3.4 Experimental and numerical investigations 
3.4.1 Experimental studies 

Despite the vast literature devoted to the experimental investigation of turbulence, very 
little is known experimentally concerning collisions of particles suspended in turbulent 
flows. 

The investigation of turbulent flows has rested for a long time on methods, such as hot¬ 
wire anemometry, which provide only information on the velocity and its spatial correlation 
function. Over the past decade, new methods have been developed, based on following 
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particles in a turbulent flows using fast-imaging [321 SO] • This has led to a wealth of new 
information on the motion of particles in a turbulent flows m- Although in principle 
feasible, detecting collisions between small particles in a well-controlled laboratory flow 
has so far not been possible. It is to be expected that this problem will be solved in a 
near future. Labelling liquid droplets with chemicals that produce an optically detectable 
reaction product when the droplets coalesce is a promising approach. 

3.4.2 Numerical investigations 

The difficulty in obtaining experimental results on collisions in turbulence makes numerical 
investigations an essential tool. The investigation of fundamental issues in turbulence rests 
on direct integration of the Navier-Stokes equation, m and (fT5]l . in a triply periodic 
domain (effectively a torus) [4111321143j . In such a configuration, the Navier-Stokes can 
be efficiently integrated using pseudo-spectral methods, based on a (truncated) Fourier 
series decompositions of the velocity field u. For efficiency purposes, the calculation of 
the nonlinear term (u ■ V)u is carried out in real space, using fast-Fourier methods to 
transform between real-space and Fourier representations. 

Once a solution of the Navier-Stokes equations is determined, the motion of parti¬ 
cles, m, can be efficiently solved m, permitting to detect collisions [45]. The available 
computer resources permit to address most aspects of particle collisions. 


4 Effect of turbulence on the collision rate 

One may think that the effect of turbulence reduces to an effective ‘rate of strain’. The 
result is more complicated, and much more interesting, as discussed in this section. In 
the following we describe various effects which become significant as the Stokes number 
is increased. In the following sections, we focus on the case of monodisperse suspensions, 
and discuss the collision rate R per particle. 

As well as influencing the settling of particles from a fluid suspension by facilitating 
collisions and aggregation, turbulent motion can have a direct effect upon the settling rate 
of heavy particles |46j . These single-particle effects are usually less significant than the 
collisional processes which are the focus of this review. 

4.1 The Saffman-Turner limit 

Particles with a sufficiently small inertia (small Stokes numbers) essentially follow the 
flow (independent of their shape or material density). In this case, the velocity gradients 
generated by turbulence strongly enhance the relative motion between particles, therefore 
enhancing the chance of collisions, by the mechanism shown in the simple example treated 
in 12.2.31 A formula derived by Saffman and Turner |8], equation ()25p below, determines 
the collision rate in the small Stokes number case, in the case where the particles are 
spherical. Two small particles of radius a advected by the flow collide if their separation 
falls below 2a. In the case of a suspension of particles uniformly distributed throughout 
the fluid with density n, the rate of collision R for a single particle is obtained, in the 
spirit of ([8]), by integrating the relative radial velocity Avr over the surface of a sphere of 
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radius 2a: 


( 21 ) 


R = -n j dS^ • lAur 

In this case the relative velocity is determined by the linearisation of the flow field, so 
that the relative motion of two particles may be described by a hyperbolic velocity field, 
such as that illustrated in figure [5l The factor of ^ in (|2T|) is required because only half 
of the sphere where the relative velocity is negative contributes to the collision rate (the 
overall flux is zero, due to incompressibility). If the particles are small compared to the 
Kolmogorov length of the flow, the relative velocity resulting from the action of the local 
velocity gradient A is given by Eq. Q. Using the explicit expression Aij = duijdxj leads 
to the following expression for the velocity gradient: 

(|Aur|) = 2a (|n • A • n|) (22) 


where n is an arbitrary unit vector on the unit sphere. Because of the isotropy of the 
velocity field, we have 

(|n-A-n|) = ( 

This leads to: 


) • (24) 

In order to estimate the expectation value of the partial derivative \dux/dx\ we assume 
that the elements of A are Gaussian distributed, and note that the variance of duxjdx 
is given by e/15z/ (see equation (fT^ l. The averaged value is equal to the inverse of the 
Kolmogorov time scale, up to numerical prefactors [8], leading to: 


R = 27r(2o)^re( 


dux 


dx 


dux 


dx 


(23) 


/Stt n(2a)^ 

a — \ - 

V 15 tk 

This estimate of the collision rate is based on the only assumption that particles are uni¬ 
formly distributed, and can be easily extended to suspensions of particles with a dispersion 
of size [8]. 



4.2 Clustering in turbulent flows 

When the effects of inertia become sufficiently large, i.e., when St is large enough, other 
effects arise, and the collision rate cannot be simply understood in terms of relative motion 
in the fluid. At finite St it is known that particles can show a pronounced clustering. We 
assume here that the local particle density remains small enough, to prevent feedback from 
the particles on the flow. 


4.2.1 Clustering and the centrifuge effect 

Particle clustering is expected to increase the collision rate. The effect of clustering is 
often ascribed to particles being expelled from vortices by centrifugal action, as proposed 
by Maxey m, based on the following argument. Integration of the equation of motion 
([I]) gives 

v{t) = I dt' u{x{t'),t') ex.p[—{t — t')/Tp] . (26) 
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(27) 


The expression (|26p reduces, in the case where Tp is small, to: 


V = u{x{t),t) 


- T, 


Du , , , , 


where D/Dt = d/dt + u . The important remark is that the velocity field v, which 
transports the particles, differs in an essential way from the flow velocity field, u: v is 
effectively compressible, with a divergence 


V 


E / dui duj \ 
\dxj dxi) 

7.1 ^ J / 


-Tptr(A2) 


(28) 


where A is the velocity gradient tensor, 
strain rate tensor (symmetric) and = 
to: 

V - V = -Tp 


Expressing A = E + fi, where = E is the 
-Cl is the vorticity tensor (antisymmetric) leads 


tr(E2) 


1 

-CJ • u 

2 


(29) 


where cj is the axial vorticity vector corresponding to the antisymmetric tensor fl, so 
that the particle flow is contracting in regions of high strain and expanding in regions 
of high vorticity. Numerical evidence of a negative correlation between particle density 
and vorticity when St ~ 1, as predicted by the ‘centrifuge’ argument has been found 
numerically |48j . We stress that depends upon making the approximation that Tp is small. 
Using (I29p provides to a method, valid in the limit of small Stokes number, to analyse 
clustering [Ml [501 EH- 


4.2.2 Clustering as a manifestation of a fractal distribution 

Figured] illustrates the type of clustering phenomena which can occur in turbulent flows. 
The clustering effect is much stronger than the simple ‘centrifuge effect’ argument suggests, 
and the positions where the greatest enhancement of density will occur are not predicted by 
that argument. In reality, the clusters have fractal properties, which should be quantified 
to determine their influence on collision rates. One way of characterising the fractal 
distribution of particles is via the correlation dimension [521153j . To define this, pick a 
particle at random, and determine the number J\f(e) of other particles within a ball of 
radius e centred on this particle. The set can be characterised by averaging AA(e). The 
quantity is related to the pair correlation function C(e) of the particle distribution: 


C(6) 


1 d(AA) 
47rne^ de 


(30) 


If the particle positions sample a fractal measure, it is expected that (AA)(e) has a power- 
law dependence: 


(W(e)) ~ nrf 



(31) 


where D 2 is the correlation dimension of the set. 
function is 


This implies that the pair correlation 
D 2-3 


C(Ar) 
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(32) 







The existence of a fractal measure is implied by very general arguments from dynamical 
systems theory |52l [53l [Ml I55j . The difference between the dimension of the attractor, 
D 2 , and the dimension of space, d = 3, is sometimes termed the dimension deficit. This 
quantity has been investigated numerically, and found to be a function of the Stokes 
number, reaching a maximum value of approximately 3 — D 2 ~ 0.7 at St ~ 0.7 [56]. The 
effect of particle clustering is argued to enhance the particle collision rate (|25p by a factor 
C(2a) (which is large compared to unity when a <C see Section 15.21 for more precise 
estimates). As well as numerical investigations, it is also possible to apply analytical 
techniques to calculate fractal dimensions (57] [58] . 

We mention here that alternative methods have been proposed to characterize cluster¬ 
ing [59] ED]. 

4.3 Caustics 

When inertial effects are more significant, particles can move independently of the fluid, 
and thus collide with a large relative velocity with other particles. As we explain here, this 
implies the existence of singularities in the phase-space of the suspended particles which 
are analogous to ‘caustics’ in optics. These caustics can lead to a pronounced increase in 
the collision rate as the Stokes number increases. It was realized [50] that when inertia is 
large enough, particles ejected from vortices acquire velocities which can be very different 
from the fluid velocity, hence run into other particles with large relative velocity Au. 
This was termed the ‘sling effect’ [50]. Later the phenomenon was described in geometric 
terms, using the known notion of caustics m, which provides a convenient framework for 
understanding and generalisation (caustic structures had been observed earlier in models 
of particle suspension El]). Another way to discuss the phenomenon has been proposed 
in [621E31E1]. 

In order to explain the role of caustics, consider a one-dimensional model, with a fluid 
velocity u{x,t), and where virtual particles, which are free to pass through each other, 
evolve with velocity v according to ([T]). Consider a cloud of particles which are initially 
on a manifold in this phase space (see figure El(a))- In the case of particles of negligible 
inertia, the velocity v of particles remains in the neighborhood of u{x,t), so the velocity 
difference between two closely separated points remains very small. If the particle inertia 
is significant, however, the velocity of particles can become multivalued (different values 
of V at the same position) because faster moving particles can overtake slower ones. This 
is illustrated in figure E](b). The points where the projection of the phase space-manifold 
onto the x-axis is singular are the caustic singularities. The caustics are focal points 
where particles with different velocities are brought together at a single point. They are 
completely analogous to optical caustics, where light is partially focussed onto a line. 
Note that caustic singularities are created in pairs, and that between them, the velocity 
is triple-valued. The existence of finite velocity differences between particles at the same 
location facilitates collisions. 

In two dimensions the caustic points are replaced by caustic lines, where the density 
of particles diverges, illustrated in figure [2] In three dimensions the caustics are sheets. 
In addition, higher-order singularities of the phase-space manifold, classified using catas¬ 
trophe theory [65], can also form. For example, pairs of caustics can originate from a cusp 
singularity. 
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Figure 6: Phase space for a one-dimensional model: (a) The particles are initially on a 
manifold with a single-valued velocity (b) If particle inertia is significant, particle positions 
may cross, because fast particles overtake slower ones. The singularities where ^ diverges 
are termed caustics. Creation of caustics is associated with the velocity becoming multi¬ 
valued. 

Caustics can increase the collision rate in two different ways. Firstly, the particle 
density becomes singular in the vicinity of the caustics at xq- If xq is the position of 
the caustic, the density has a generic divergence of the form n(x) ~ (x — xo)~^^^ as 
the caustic is approached from one side (here x is any generic coordinate). Secondly, as 
already explained, the multivalued character of velocity can greatly enhance the collision 
rate. The same mechanism has also been discussed in terms of ejection from vortices [50], 
and termed sling effect. The relations between caustics and collisions of inertial particles 
were proposed in [HI |66| . 

The nature of the singularity at caustics can be illuminated by considering the gradient 
of the velocity for nearby particles. Assuming a continuum fluid description for the flow 
of particles, and differentiating Eq. ([T|) leads to the following evolution equation for the 
tensor cr, defined by aij = diVj\ 


d(T 9 1 / . \ \ 

—+cj2 = —(A-Cj) . (33) 

Oju Tp 

where A is the velocity gradient tensor. Equation (1331) . whose nonlinearity is of the Burgers 
type, is a convenient starting point to analyze caustic formation numerically m or even 
experimentally [68j . 

To estimate the contribution of the caustic mechanism to the collision rate, Eq. ([8]). 
we approximate the typical velocity difference associated with the sling effect as v/tk, 
up to a dimensionless function of St and the Reynolds number of the flow. Re: 

fl = ^^F(St,Re). (34) 

tk 

The function F(St, Re) must vanish in the limit as St —)• 0, because the caustic mechanism 
is absent in the advective limit. In the Appendix it will be argued that St ^ 0 is a singular 
limit, and that 

F(St, Re) (X exp(—S'/St) (35) 
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provides a qualitatively plausible representation, in agreement with numerical estimates ISZl 
[69]. 

It will be demonstrated later, see Section |5l that the caustic or sling effect provides 
the dominant contribution to the collision rate, as soon as the inertia becomes significant 
(for St > 0.5). 

4.4 Random uncorrelated motion 

In the limiting case where the turbulence intensity is very high, an alternative approach 
to understanding the effect of increasing the turbulence intensity was initiated by Abra- 
hamson m- When the damping timescale Tp is large, particles travel for long distances 
through the fluid, without being influenced by the smaller scale turbulent eddies. Abra- 
hamson argued out that this implies that the particles which reach any one point will have 
approximately isotropic random velocities, so that the gas-kinetic approach described in 
section (2.2.11 can be used to model the motion of the suspended particles. This point has 
also been emphasised in |63j . where this regime was termed ‘random uncorrelated motion’. 

In a turbulent flow there is a largest timescale, rx, (the integral timescale) which is 
characteristic of the driving motion which creates the turbulence. Whenever Tp <C rx, 
the particles will be advected with the largest eddies, and not be sensitive to eddies with 
a timescale smaller than Tp. The velocity difference (|Au|) in Eq. (Hjcan be estimated as 
the typical relative velocity of particles at any given point. This is comparable with the 
relative velocity of the particles and the fluid. We now describe a simple argument which 
can be used to surmise an expression for the mean magnitude of the relative velocity, 
(|Au|). 

We argue, in the spirit of Kolmogorov phenomenological description, that in the inertial 
range the only parameter which enters the description of the flow is the rate of dissipation 
per unit mass, e, together with the time scale characteristic of the particle motion, Tp. 
Dimensional arguments then lead to the following expression for the relative velocity: 

{\Av\)=IC^ (36) 

where /C is a dimensionless constant (this approach was developed in [71]). This leads to 
a formula for the rate of collision in a highly turbulent flow: 

R^K^VSt (37) 

''“k 

The constant K, which is potentially important for astrophysical applications, will be 
discussed in Section 15.21 

5 Synthesis: a unified theory for the collision rate 

5.1 Modelling the collision rate 

We have discussed several mechanisms for the role of turbulence in facilitating collisions 
between suspended particles. The shear component of turbulence causes collisions between 
particles moving with the fluid, at the rate which estimated by Saffman and Turner. In 
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addition, the inertial effects, parametrized by the Stokes number, induce two distinct ef¬ 
fects. Firstly, there is the clustering effect known as preferential concentration. Secondly, 
particles located at the same spatial location may have very different velocities, a conse¬ 
quence of the caustic folds in phase space. Finally, in the limiting case when the inertia is 
very large, it was argued that the motions of the particles become completely uncorrelated 
from one another, and a simple asymptotic form for the collision rate was proposed. In 
this section, we show how to combine these competing mechanisms. The main result of 
the section is a derivation of a simple expression for the collision rate. It involves some 
parameters, which must be determined empirically by comparison with simulations. We 
also consider the evidence supporting our expression coming from numerical simulation as 
reported in m 

The central idea is that the collision rate can be resolved into two components. Some 
collisions arise from particles which follow similar trajectories for an extended period, and 
which eventually come into contact because of shearing motion in the flow. At low Stokes 
number, where particles are exactly following by the flow, this is the only mechanism. The 
collision rate due to this mechanism is denoted i?adv j and it is controlled by the local shear 
rate, as well as by the local concentration around particles. Collisions between particles 
not following fluid path lines, occurring when caustics start to form in the phase-space of 
the suspended particles, give rise to a very different contribution, denoted here by i?caust- 
The two contributions i?adv and iJcaust differ in an essential way by their dependence on 
the siza a of particles. We argue that these mechanisms operate independently and that 
their contributions are additive, so that 


R — I?adv + R 


caust • 


(38) 


This decomposition was proposed by analysing theoretical models [66l Ea] and simula¬ 
tions of simplified numerical models [69]. There is no hard criterion which distinguishes 
between an advective and a caustic-mediated collision, so (f38ll must be regarded as an 
approximation. In the following, we derive the form of the terms iiadv and i^caust in (jMj)- 
In the limit St <C 1 the collision rate is determined by shearing motion, without any 
inertial effect, and the collision rate is given by equations (ESI) E- 

The preferential concentration causes clustering of particles with finite values of St. 
The density of particles at a distance r from a given test particle is nC{r), where C{r) is 
a radial correlation function. An important remark is that the inertial effect that leads 
to clustering does not enhance signihcantly the relative velocity of particles, contrary 
to caustics. In fact, particles with large velocity differences do not stay together, thus 
not contributing to preferential concentration. As a result, the enhancement of the local 
concentration around particles contributes only the i?adv term, which becomes: 


.^adv 


87rn(2a)3 



tk 


C(2o) 


(39) 


When inertial particles converge to a fractal measure, as we have argued, see Section. [4.211 . 
the function C{r) behaves as a function of r as a power-law: C{r) oc 

In the caustic-dominated case we combine expressions (IM]) . ([35]) and (l37|) to obtain an 
expression for the rate of collisions due to caustics: 

2 

^caust = K^V^exp(-5/St) (40) 

Tk 
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where K and S are two constants to be determined by fitting the collision rate against 
numerical simulations: K determines the asymptotic collision rate at very large Stokes 
numbers, and S determines the crossover point where the caustic mechanism becomes 
significant. The non-analytic term is consistent with theoretical expectations [66] and 
numerical studies m- 

Taken together, equations (l38|) - (H0]) define the theory for the collision rate. We note 
that the function F{St, Re) defined in equation (|.S4jl is known only through its asymptotic 
behaviour in certain regimes, so further information is necessary to make progress. The 
two terms in (|38p differ in an essential way through their dependence on the size of the 
particles a, at a fixed value of St and Re. Namely, Rcaust varies as a^, wheres i?adv varies 
as {D 2 > 2). This difference can be used to separate the different contributions 
numerically. We present in the following section numerical results, which indicates that 
the decomposition (I38p is indeed a very effective tool for the analysis of collision rates. 
We note that Zaichik and co-workers [Ill [75] proposed models of the rate of collision 
in turbulent flows which use similar physical principles. Their final expressions for the 
collision rate involve a much larger number of parameters than equations (|38p - (|4n|) . 

As well as the collision rate, the probability distribution of relative velocities of colliding 
particles is also of interest, especially in the context of understanding planet formation, 
where collisions may be sufficiently energetic to cause fragmentation [761177]. 

5.2 Numerical evidence on collision rates 

A substantial amount of literature has been devoted to the numerical investigation of the 
collision rates in both direct numerical simulations (DNS) of the Navier-Stokes equations, 
and model flows, but many of these studies pre-date some of the theoretical insights 
contained in equations (j38[) - ()4n[) above. The collision rate does show a marked increase as 
the effects inertia are increased, and this is usually ascribed to the effects of ‘preferential 
concentration’, that is the clustering effect, but it has been argued that effects of caustics 
may also be significant [30]. Here we summarise some recent numerical studies which 
separate out the two contributions in equation ((38P by studying the dependence of the 
collision rate upon particle size, keeping the Stokes number fixed. This is achieved by 
varying the ratio of the particle density to the fluid density, Pp/pi- Modifying the ratio 
Pp/Pi at fixed value of the Stokes number is achieved by varying in the collision detection 
algorithm the radius of the particles, a, according to (d,® (so that a oc {pp/Pi) ^^^). 

The collision rate, i?, is determined numerically by recording among a set of trajecto¬ 
ries, all instances in which the separation radius decreases past 2a. In the case of collisions 
where particles stick or coalesce on contact, we should only count the first contact col¬ 
lisions. This effect should be accounted for by introducing a factor / < 1 in (I38p . The 
coefficient is no smaller than ~ 0.85 when St is very small, and decreases when St in¬ 
creases m- Here, we do not distinguish between single and multiple collisions. The 
collision rate, R, shown in Fig. [TKa), is normalized by n(2a)^/rK and plotted as a function 
of St. The Saffman-Turner prediction, (1251) . implies that in the limit St —)• 0, the quantity 
RrK/(n(2a)^) should become independent of the ratio Pp/pi- Our own numerical results 
m, in agreement with other estimates |30( [3T] , are only consistent with this prediction for 
small values of St. Fig. db) shows that Rtk/{ na?rf) as a function of the Stokes number, 
does not depend much on Pp/pi for values of St larger than > 0.3. This scaling is con- 
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Figure 7: The collision rate R as a function of the Stokes number St, for three different 
values of the ratios of density Pp/pi = 250, 10^ and 4 x 10^. The collision rate R is 
normalized by n(2o)^/TK (a), and n(2a)^7//rK (b). The horizontal dashed line in (a) 
corresponds to the Saffman-Turner prediction. 


sistent with the sling/caustics collision mechanism, described by equation We note 

that F{St, Re) deduced from Fig. [Tl^b) does not fit the asymptotic form F(St, oo) = R-v/St 
for large values of St. We ascribe this to the limited Reynolds number of our numerical 
simulations. 

Figure [8] shows the effect of clustering in the same simulations: the effect of clustering 
on the function C{2a) reaches a maximum at a value St ~ 0.7 of the Stokes number. The 
collision rate grows much larger at higher values of the Stokes number, which provides 
further indication that caustic effects dominate the collision rate. 

Figure [9] shows the relative importance of the advective and caustic terms in these 
simulations. The caustic contribution becomes dominant for Stokes numbers greater than 
« 0.75. 

At very large Stokes numbers the collision rate is expected to be given by (13711 . which is 
a potentially important prediction for astrophysical applications. To estimate the unknown 
constant K in (j37ll . one needs simulations at high enough St, i.e., high enough Tp, but with 
the constraint that the value of \R\, such that Tp = (|i?p/e)^/^, is in the inertial range. 
Extrapolation of numerical data at moderate Reynolds numbers from different sources m 
leads to a value 77 Ri 50. 

An alternative decomposition of the collision rate, originally proposed by m , expresses 
the collision rate R as a product in which the term C{2a), which describes the local 
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Figure 8: The function C(2a) that measures preferential concentration, for three values 
of the density ratio, Pi/pu as indicated in the hgure. The preferential concentration does 
not play a significant role for St > 5 



Figure 9: The ratio of the contribution to the collision rate due to preferential concen¬ 
tration -Radv) defined by equation (l3^ and of the total collision rate, R. At fixed value 
of Pp/Pi, The contribution of i?adv to the total collision rate decreases when the Stokes 
number increases, or when the ratio Pp/Pi increases. 
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concentration enhancement around a particle, appears as an overall factor: 

R = 27r(2a)^nC(2a)(| AuDeff (41) 

This representation, which is exact for a suitable definition of (|Au|)efr, suggests that 
the preferential concentration and sling effects act together to enhance the collision rate. 
Fig. [71(b) demonstrates that if this parametrisation of the collision rate is used, then the 
dependence of C(2a) upon Pp/pf, see Fig. [8]must be cancelled (for St > 0.5) by a reciprocal 
dependence of the collision velocity, (| AuDeg. In fact, previous measurements [HOllHT] of the 
dependence of C{r) and of the average velocity difference as a function of r suggest power 
law dependences, the exponents being such that the product C(2a)(|Au|)eff is essentially 
constant for St > 0.5. Equations (l38]l . (13^ and (140l) provide a very natural explanation 
of this cancellation. We remark that the power-law dependence of the collision velocity 
has been explained in a random flow model 182j. and used to justify Eq. (I38p for that 
system. 

We therefore conclude that the decomposition (I38p . which rests on a physically well- 
motivated analysis, and which is supported by the analysis of simplified theoretical models, 
provides a consistent description of the available numerical data. One of the main lessons 
from the analysis of the dependence of the collision rate on St and on the ratio Pp/pf, 
is that the sling/caustic effect provides the dominant mechanism for the dramatically 
enhanced collision rate of particles whose Stokes number exceeds > 0.75, in the cases of 
water droplets. 


6 From collisions to aggregation 

Up to this point we have considered expressions for the rate of geometrical collisions. In 
order to determine the aggregation of particles, we need to consider whether the collisions 
result in particles combining, and how the population of particle sizes evolves. 


6.1 Collision efficiency 

Droplets which undergo a geometrical collision, i.e. whose separation falls below oi -|- 02 
when their motion is predicted by equation ([1]), might not coalesce, because the streamlines 
of small droplets curve around larger ones. In fact, if the Navier-Stokes equations were a 
complete description, droplets would never collide, because of the presence of a lubricating 
film of fluid between them. Other mechanisms allow the particles to actually coalesce [83]. 

The aggregation of particles is characterised by a collision efficiency e, defined as 
the ratio of the observed rate of coalescence to the rate of geometrical collisions (where 
the separation of the centres predicted by Eq. © falls below oi + 02 )- The coalescence 
efficiencies are hard to determine, because the the complex physics determining breakdown 
of the lubricating layer. It is widely accepted that they are low for typical cloud droplets 
|32l|2|. If the larger droplet has radius below 20^m, it is believed that e < 0.1, and that 
for radius 10/xm, e < 0.03 |2]. Eor droplets of size a = 50/rm colliding with droplets of 
size a = 10 //m, however, the efficiencies are expected to be close to unity [32l[2|. 
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6.2 Models of aggregation 

The theory developed by Smoluchowski [12] to describe the coagulation in colloids seems 
to provide a general framework to discuss a vast class of problem, including coalescence 
of water diiHii or of dust grains. The approach of m is naturally formulated in terms 
of mass conservation, so we use here N{m, t) the number of particles of mass m per unit 
volume at instant t. The Smoluchowski equation is: 

dN If”* - - 

= 7 T / dm' K{m — m',m')N{m',t)N{m — m',t) (42) 

at 2 Jq 

poo 

— N{m,t) / dm'K{m,m')N{m',t) . 

Jo 

The Smoluchowski description is a mean-field model, in that it assumes that the particle 
density is spatially uniform. To treat aggregation processes, the Smoluchowski approach 
is generally preferred to the full master equation, which is completely intractable. 

Depending on the structure of the collision kernel K(m.m'), the Smoluchowski equation 
can predict one of two different types of particle growth. In ripening processes, the particles 
grow but they remain of comparable sizes. Alternatively, if the collision kernel K{m,m') 
increases sufficiently rapidly as m, m' oo, the growth process is unstable: a small 
fraction of particles undergoes runaway growth in a finite time. This means that, in the 
Smoluchowski model, all of the smaller particles end up absorbed into one massive cluster 
in a finite time. This phenomenon, known as gelation, has been investigated intensively: 
see m for a review. In the notation of Eq. (|42p . gelation is expected to happen when the 
collision rate K{m,m') grows faster with m and m' than linearly [85]. This is of primary 
concern in our case, since for large values of the mass, or equivalently, of the radius, the 
dominant physical effect in the collision rate is due to particle settling, see Eq. (1121) . In this 
case, the kernel K grows as K{M, m) oc when M is large, so the mean-field equation 
([42]) leads to gelation. Here we discuss the consequences of the gelation transition in the 
Smoluchowski model. 

The requirement of a very large number of collisions, of the order of 10®, to form a 
raindrop of typical size 1 mm out of much smaller microscopic droplets (of typical size 
10 /rm) makes gelation an appealing feature of the Smoluchowski equations. 

There is, however, a potential problem with using the Smoluchowski equation to model 
rainfall. The gelling transition occurs after a critical time tc- Investigations of the gelation 
transition for the Smoluchowski equation for homogeneous kernels of the form 

K{mi,m 2 ) = const, m'^m'i (43) 

have shown rigorously |86l 187] that there is a gelation transition when + v > 1, and that 
tc = 0 (termed instantaneous gelation) whenever v > 1 or /i > 1. The kernel for liquid 
droplets settling under gravity is asymptotic (when mi ^ m 2 ) to a homogeneous kernel 
with V = 4/3, so the issue of instantaneous gelation is relevant. 

These results suggest that the Smoluchowski equations, which predict an instantaneous 
gelation, must to be used with caution. 

We should consider the reason why the Smoluchowski equation has this pathological 
behaviour when u > 1. Two observations are necessary to explain the origin of the 
zero-time singularity. Eirst, note that a cluster of size M will grow to infinite size in a 
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time which should be a decreasing function of M (because the rate of collisions 

is assumed to increase with the cluster size). Furthermore, under the conditions where 
instantaneous gelation occurs, this time tc{M) approaches zero as M —>• oo. The second 
observation is that, in an infinite system, a cluster of size M arises in a finite time to 
somewhere in the system, no matter how small we choose to- can find a cluster of size 
M somewhere in an infinite system, no matter how small we make to or how large we make 
M. According to the Smoluchowski equation, the model thus undergoes gelation at a time 
t = to + tc{M) where both to and tc can be made arbitrarily small. This argument applies 
to a system with an infinite number of particles: see [88] for a discussion of finite-size 
effects. 

We note that the prediction of a runaway growth may need to be revised, because 
when the droplets become very large, the terminal velocities are proportional to so 
the collision kernel grows as m®/®, which does not yield a gelation transition. In addition, 
rain droplets may undergo fragmentation when they become sufficiently large [89]. 

This pathology of zero-time gelation results from the use of a mean-field description, 
which ignores any spatial structure of the particle density N(m,t). Under the mean-field 
approximation, the local depletion of particles due to the formation of a large cluster 
is not taken into account, so an unbounded growth can proceed. We conclude that the 
Smoluchowski equation must be used caution, and that a proper description of the runaway 
growth therefore requires more realistic assumptions about the spatial distribution. 


7 Application to rain initiation 

The problem we are focusing on here concerns the generation of rain drops in ‘warm’ (ice- 
free) cumulus clouds [32l [2l [90[ [3] . In ice-bearing clouds the Wegener-Bergeron-Findeisen 
can result in rapid growth of ice crystals by condensation, which makes it easier to ex¬ 
plain precipitation [32]. When the air becomes supersaturated, water droplets condense 
rapidly onto aerosol nuclei. They grow up to a typical radius 10 fim before the supply 
of water vapour is exhausted. The later stages of the growth of a raindrop (from radius 
of approximately 50 jim upwards) involve the falling drop growing rapidly by coalescence 
with microscopic droplets lying in its path, and are easier to understand. The challenge 
is to explain the growth of droplets from size 10 fim to approximately 50 /rm, and in par¬ 
ticular, to identify the physical process leading to the rapid onset of rain showers, which 
can develop in less than half an hour. 

To illustrate the discussion, consider the following representative values for a convecting 
cumulus cloud which could produce precipitation The typical droplet radius is 

a = 10/xm, the number density of droplets is re = 4 x 10® m“^, and the cloud depth is 
L = 10® m. The typical vertical velocity of air inside the cloud has magnitude 2ms“^, so 
that the eddy turnover time may be taken to be tl = 10® s. An estimate for the rate of 
dissipation is e ~ 1?/t\ = 10“®m^s“®, which gives an estimate of the Kolmogorov time 
tk ~ 10“® s. Rain falls as droplets of size approximately 1 mm. 

What are the processes possibly leading to the growth of very small droplets with an 
initial size of typically 10 //m? Growth can proceed by collisions, and one possibility is 
that droplets settling at different rates lead to collisions. This mechanism is only possible 
in the presence of size dispersion among droplets. Specifically, inserting values for air and 
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water at 5°C into equation (1111) gives k ~ 1.4 x 10^ Then, using the expression of 

the collision rate (fT^ to estimate the coalescence rate of a droplet of radius a + Aa (with 
a = 10 fim, Aa ~ 2.5/im) falling through a gas of smaller droplets (of radius a = 10/rm), 
and with a collision efficiency e ~ 0.03, gives a rate of coalescence R ~ 10“‘^s“^. We 
conclude that the rate of coalescence of typical sized water droplets induced by differential 
settling is very small, as long as the dispersion between particles is small. 

As the air in a convecting (cumulus) cloud is turbulent, it may be expected that the 
coalescence rate between droplets may be greatly facilitated by turbulence [8]. Using the 
expression for the collision rate between small droplets following the flow, (I25p . and the 
parameters of the cloud model, gives i?turb ~ 2 x 10“® s“^, which is negligible. The effects 
of turbulence are dramatically increased when the effects of droplet inertia are significant. 
Inertial effects are measured by the Stokes number, St = Tp/rx. The collision rate is 
greatly enhanced by effects due to caustics for St > 0.3 m- With the values of cloud 
model given above, however, the estimated Stokes number is St ~ 10“^, which is to small 
to lead to any significant enhancement. 

The small rate of collisional coalescence of droplets is clearly a problem. It should 
be kept in mind that considerable variation exists among various clouds, and even within 
a given cloud, the flow is expected to be very inhomogeneous. The numbers used in 
the estimates above could therefore vary significantly, resulting in a large increase of the 
coagulation rates, compared to the values given here, at least in parts of the clouds. One 
may also remark that only a very small proportion of the microscopic droplets needs to be 
converted into raindrops. Consider the rate at which droplets actually form. Rainfall at a 
rate of 3.6 mm hr“^ = I0“® m s~^ is considered as ‘moderate’ [32] . If the raindrops have size 
a ~ 1 mm, the number of drops falling per second and per square meter is approximately 
250. Given the assumed cloud depth of L = 10^ m, the volumetric rate of production 
of raindrops is approximately 0.25m“^s“^. If the microscopic droplets have density n = 
4 X 10® m~^, then the rate of conversion of each microscopic droplet into a ‘collector’ 
droplet undergoing runaway growth is approximately 6 x 10“^® s“^. Alternatively, during 
a five minute shower, the probability that a water droplet starts growing, and accumulates 
during its subsequent fall enough droplets to become a rain droplet is small, approximately 
2 X 10“^. The problem of rain initiation is, therefore, concerned with the frequency of 
very rare events m- 

Despite the fact that the required conversion probability is very small (of order 10“^), 
growth of droplets is too slow by a collisional mechanism. On growing from 10 fim to 
50/nn, the volume of a droplet increases by a factor of 125, that is, there are or order 
100 collision events. It was argued above that the rate for the first collision events is 
small, R ~ 10“^ s“^. It is not obvious whether the rarity of the event is sufficient to 
compensate for the low rate of multiple collisions. Large deviation theory (reviewed in 
I92|) is the appropriate tool for analysing this problem. A large deviation analysis shows 
that a sufficient number of droplets can undergo runaway growth in small fraction of the 
mean time to the first collision |93j . 

After a droplet has grown to a size much larger than the typical droplet size, it falls 
rapidly and collects other droplets in its path. When a > 50 ^m we assume that the 
collision efficiency is approximately unity. Consider a droplet of size oi falling through 
a ‘gas’ of small droplets, which can be characterised by the liquid volume fraction Q = 
ATTn{a^)/3. The large droplet falls with velocity v = Ka\ and grows in volume at a rate 
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TrafQv, so that 


(44) 


doi kQoi 
~dt ~ 4 ■ 

Solving this equation shows that the droplet radius diverges in the time 

4 

kQui 


(45) 


Equation (I45p predicts that the time before runaway increases rapidly as the droplet size 
decreases. For the parameters introduced to describe a cloud, a droplet of size oi = 50 /rm 
requires Texp ~ 2 x 10^ s to undergo explosive growth. 

We conclude that, although runaway growth can proceed when a droplet reaches 50 /xm, 
it is difficult to understand how droplets can reach this size in the short time that it takes 
for a rain shower to develop. While turbulence mediated collisions may be an important 
ingredient to explain the runaway growth of rain drops in warm cumulus clouds, a complete 
picture is likely to involve other effects, possibly including those triggered by non-collisional 
growth processes [94]. Such effects need to be discussed more systematically in the future. 


8 Application to planet formation 

The other major area for applications of turbulent collision processes is in understanding 
planet formation. Here we can only give a brief introduction to a complex and rapidly 
developing research field. We start by summarising the standard model, which is reviewed 
in m, before discussing the unresolved difficulties. 

When a star forms by gravitational collapse, conservation of angular momentum pre¬ 
vents all of the cloud of gas from falling into the star, and the residual material rapidly 
forms a disc-like structure (which has the least internal motion consistent with the con¬ 
servation of angular momentum). It is assumed that planet formation occurs in such 
circumstellar discs surrounding young stars [7|, which implies that planets would be cre¬ 
ated with near circular orbits in the plane of the disc. This model is consistent with the 
structure of the solar system, where the planets lie in approximately circular and coplanar 
orbits, close to the equatorial plane of the Sun. The structure of these circumstallar discs 
is described by a model introduced by Shakura and Sunyaev [96] . 

The interstellar medium from which planets form is thought to contain sub-microscopic 
dust. The grains have a broad distribution of sizes, but 10“^ m is a typical value. The 
grains may be ice particles or minerals, predominantly silicon or carbon based, originating 
from nuclear reactions in stars. The proportion of ‘heavier’ elements (that is, elements 
other than H and He) in star forming regions is of the order of one percent by mass. 
There is also a consensus that the dust grains play an important role in planet formation. 
It is assumed that dust grains adhere on contact, making ever larger structures. In the 
simplest version of this model, these objects continue to accumulate material until they 
become kilometre-sized ‘planetesimals’, and gravitational forces between the planetesimals 
take over. A variation on this model suggests that ‘boulders’ settle to form a layer at the 
mid-plane of the circumstellar disc. When the density of material at the mid-plane reaches 
a critical level, this triggers a gravitational collapse [97] . 
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Collisions between grains to produce larger structures play an essential role in these 
routes to planet formation, and the aim is to understand whether a quantitative descrip¬ 
tion of these processes is viable. In addition to the difficulty in estimating the physical 
parameters inside a circumstellar with current observational techniques, some serious the¬ 
oretical difficulties with these models are hard to circumvent [98]. 

In the absence of any turbulent transport mechanism, the rate of collision of micro¬ 
scopic grains is very low. It seems likely, but it is not certain, that the gas in a circum¬ 
stellar disc is in turbulent motion: the Reynolds numbers associated with the motion are 
extremely high, and the effects such as convection due to frictional heating will combine 
with the rotational motion. Turbulence certainly increases the rate of collisions, but it 
also brings its own problems |99] . The dust particles will adhere to each other due to van 
der Waals forces and to a certain extent electrostatic forces. Because these binding forces 
are weak, aggregates of dust particles are very easily fragmented by collisions. Estimates 
of the relative velocity of particles colliding in a turbulent environment such as equation 
(I36p indicate that the collision speed increases with the size of the particles jH [7T], so 
that there may be a maximum value for the size of a dust cluster which can be formed by 
aggregation in a turbulent environment. Recent estimates indicate that the maximum size 
that can be reached by clusters of dust particles appears to be very small for reasonable 
values of the parameters in a model for the protoplanetary accretion disc. It is useful to 
give some estimates for conditions inside the circumstellar disc according to the Shakura- 
Sunyaev model. These have a power-dependence upon distance from the star, but at 
lAU = 1.5 X 10^^ m (i.e., the Earth-Sun distance) the gas density is 3 x 10 ^ kgm the 
dissipation rate is 10“^m^s“^, the gas mean-free-path is 4 x 10“^ m, the speed of sound 
is 700 ms“^ [98|. The low density of the gas implies that the Epstein formula ([3|) is appli¬ 
cable, and the particle relaxation time Tp is very large. Thus, equation (f36l) implies that 
the relative velocity of colliding particles is quite large. For example 30 cm size objects are 
predicted to collide with velocities of approximately lOrns”^ |98] . It seems improbable 
that balls of dust particles would survive collision at these speeds. These considerations 
suggest that there are severe theoretical difficulties with models based upon aggregation 
of dust particles. 

Other outstanding theoretical difficulties are related to the model itself, and are inde¬ 
pendent of whether turbulence plays a role in collisional aggregation. One of the difficulties 
concerns the fact that the gas in an accretion disc is partially supported by its pressure, 
so that its orbital velocity is slightly lower than the Keplerian value in the quasi-static 
state. A ‘rock’ (more accurately, an aggregate of dust grains and possibly ices) which is 
entrained with the gas is not supported by the pressure and consequently slowly spirals in 
towards the star |in0[ IlDlj . This effect is most pronounced for rocks with size comparable 
to the mean free path of the gas (typically 1 cm to 1 m), and the timescale for spiralling in 
is of the order of 300 yr starting from an orbit at lAU. Even under the most favourable 
assumptions about growth rates by aggregation, it is difficult to see how aggregates of 
dust particles can grow sufficiently rapidly to avoid spiralling in. A variety of compli¬ 
cated theories have been proposed to try to resolve this difficulty. As of now a ‘streaming 
instability’ |102l 1103] is regarded as a promising theoretical approach. 

More recently, the detection of a large number of extra-solar planets |104j has brought 
new challenges. These discoveries yielded many surprises, some of which are hard to 
reconcile with the standard model for planet formation. Significant numbers of these 
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exoplanets have large orbital eccentricity. Various models have been proposed to account 
for this |105j . The most plausible of these is a slow-acting three-body instability resulting 
in a drift of orbital parameters, leading to a near-collision between two planets. This could 
cause escape of one planet and scattering of the other to an eccentric (and probably non- 
equatorial) orbit [106] . It is as yet not clear whether the large proportion of exoplanets 
with eccentric orbits can be explained by this model. The model would suggest that 
large planets are less likely to be scattered into highly eccentric orbits. There seems, 
however, to be a positive correlation between eccentricity and mass, which lends support 
to alternative mechanisms of planet formation [107] . There are also numerous examples 
of planets where the orbital plane is at a very large angle of inclination to the equator of 
the star, or where the orbit of the star is in the opposite direction to the spin direction 
of the star m- These observations cast a doubt on one of the central hypothesis of 
the standard planet formation model [7], which imply that planets are formed in circular 
orbits in the circumstellar disc. 

As for rainfall, understanding the role of turbulence in mediating collisions has not 
resolved the difficulties in explaining planet formation. If the dust accretion model will 
ultimately be shown to be correct, collisions in turbulent flows will play an important role. 
It is possible, however, that alternative theories will be required [in7( 1108] . 


9 Perspectives for future work 

The recent developments in the theory of collision between particles in turbulent suspen¬ 
sions presented here correspond to quantitative progress in describing a complex physical 
phenomenon. The collision kernel for tracer particles which are exactly following the 
turbulent flow had been understood for many years. In contrast, the problem collisions 
between inertial particles, which is much more realistic for numerous applications, involves 
subtle physical effects, which can be traced back to the fact that the trajectories deviate 
from the fluid ones. The first effect, preferential concentration, i.e. the tendency for par¬ 
ticles to be unevenly distributed in the flow, had been noticed in many experiments and 
simulations. It can be qualitatively understood in terms of fractal attractors. The notion 
of slings or caustics, leading to strong velocity difference between colliding particles, has 
emerged from a combination of theoretical and numerical work which we have reviewed 
in this paper. Over a vast range of particle inertia, the latter effect is the dominant one 
to determine the large collision rate between particles. While the understanding of these 
phenomena was largely motivated by a problem of cloud physics, it is very likely that the 
knowledge developed here will find other applications in other helds of science. 

While the knowledge of collision kernel is a crucial first step to understand coagulation 
of particles in a turbulent suspensions, an accurate description of the actual formation and 
growth of clusters rests on the solution of kinetic models. The usual approach is a ‘mean- 
field’ description originally proposed by Smoluchowski. Given the functional dependence 
of the collision kernel in problems inspired from cloud physics, this model predicts a 
singular behavior at zero time, which points to difficulties in the use of the model. The 
runaway growth of a small minority of drops is however expected to be an important 
feature in cloud physics. Identifying the mechanisms which trigger a runaway growth, and 
quantifying how frequently this happens are significant areas for further development. 
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The knowledge gained from the description of the collision kernel in turbulent suspen¬ 
sions has not yet led to an unambiguous understanding of the formation of large clusters 
in the important physical contexts of rainfall and planet formation. This situation calls 
for future work. Future theoretical developments may be necessary, in particular in the 
description of the aggregation process and understanding the significance of the runaway 
growth. However it is also likely that new experimental work to describe cloud droplet 
aggregation [109], and new observational discoveries on exoplanets [104) will provide the 
missing insights which are required to complete the picture. 
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10 Appendix A: One-dimensional model for clustering and 
caustics 

In view of the role played by clustering effects and caustic formation, it is desirable to see 
both of these effects in operation in a simplified, analytically tractable model for particle 
motion in a random flow. This appendix describes an exactly solvable one-dimensional 
model [110] which explains why the particle distribution is a fractal [58] , and why equation 
(l35l) is a good model for the rate of formation of caustics. While the model is only 
exactly solvable in one dimension, the qualitative predictions are easily extended to higher 
dimensions. 

10.1 A one-dimensional model 

An incompressible flow in one dimension has no spatial variation, and cannot generate 
either clustering or caustics. We will consider a compressible one-dimensional model: 

X = V ; V = 'y[u{x,t) — u] (46) 

where 7 is the inverse of the particle characteristic time, and u{x,t) is a random gaussian 
velocity field, delta-correlated in time, with smooth spatial correlations: 

{u{x, t)) = 0 ; {u{x, t)u{x', t')) = C{x — x')5{t — t') 
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10.2 Preferential concentration and caustics: a general argument 


To understand both clustering and caustics, we write the equation of motion for small 
separations in positions, 6x, and velocity 6v obtained by linearizing Ea. (j46p : 

5x = 6v ; 6v = —'ydv + f{t)Sx (47) 

where f{t) = Now consider the change of variable from {5x,6v) to {Y,Z), 

defined by T = In 5x and Z = |^. In terms of these variables, the linearised equations of 
motion (j47p become: 

Y = Z (48) 

Z = -7Z + Z2 + /(t) (49) 


The variable Z obeys a stochastic differential equation, independent of Y, with / acting 
as a forcing term. The variable Y is then determined by a very simple evolution equation, 
which can interpreted as a generalised random walk. The evolution of Y is then charac¬ 
terised by a drift velocity and a diffusion coefficient. The equation describing the evolution 
of Y is invariant under the transformation Tg :Y —>• T -|- <5, so the probability density of Y 
must also be an eigenfunction of Tg. This implies that a steady-state probability density 
for Y must be of the form P(Y) = 74exp(ay), which in terms of the original variable 6x, 
leads to: 

P{5x) = 5x°'~^ (50) 

The linearized equation of motion (I48I49|) is only valid when 5x is very small, corresponding 
to sufficiently negative values of Y. The solution (1501) corresponds to a normalisable 
distribution only if a > 0. 

Equation ([50P predicts that the expected number of particles in a ball of radius e 
surrounding a test particle is {Af{e)) ~ where D 2 is the correlation dimension [58]. 
The corresponding probability density is P{Sx) ~ so using (l50P leads to the 

conclusion that D 2 = a. This simple argument, which is equally valid in higher dimensions, 
establishes why particles cluster onto a fractal set. 

Crucial to the statistical properties of 5x is the Lyapunov exponent. A, defined by: 


A = lim -In 

t—>-oo t 



(51) 


where 5x{t) is the infinitesimal separation of two trajectories. It follows that 


A = (Z) 


(52) 


The variable Z also has a clear interpretation in relation to the existence of caustics. 
Caustic singularities correspond to points where 6x = 0, so that at caustics |Z| —)• -|-oo. 
The nonlinear term in (1491) for Z leads to a finite time singularity, with Z —)• -|-oo. The 
trajectory returns with Z = —00 after the singularity. 

In summary, the equation of evolution (j49l) for Z{t) provides important information, 
concerning both preferential concentration (effectively the exponent a in ()50p i. as well as 
caustics formation, which is the rate of formation of singularities of Z. 
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10.3 Solution via a Fokker-Planck equation 

The fractal clustering and rate of caustic formation can both be analysed for the one¬ 
dimensional model (j46p . with the velocity field which is white-noise in time, ()47p . Here, 
we concentrate upon the calculation of the rate of caustic formation, in order to show 
that this has a singular behaviour, analogous to (1351) . in the limit where inertial effects are 
negligible. The Lyapunov exponent, as well as the correlation dimension can be accurately 
determined for this model |110j . 

In our model, the stochastic term / in (I49p is also white-noise, with zero mean, and a 
variance {f{t)f{t')) = 2D6{t — t'), the diffusion coefficient D being simply deduced from 
the correlation tensor of u: 

O = -5^(0) . (53) 

The probability density for Z is determined by a Fokker-Planck equation |110] 

This is in the form of a continuity equation. Steady-state solution in one dimension 
correspond to a uniform flux = J, which determines the rate of escape of Z{t) to infinity, 
hence to the formation of caustics. The search of a solution is facilitated by introducing a 
dimensionless variable z, a dimensionless parameter e, and a potential (p{x,£) 



In terms of the dimensionless variables, the distribution of the scaled variable z is m 

P{z) = [ exp[(/)( 2 :',e)] dz' (56) 

^7 7-00 


Imposing the normalization condition for the probability density leads to an expression for 
the rate of caustic formation, J. In the limit as e —0, the integral over x' is approximately 
\/^exp(I/6e^) when is in the interval [—1/e, l/2e]. Integrating over Z to normalise the 
distribution yields the approximation 


J = 


SL 

2tt 


exp 


1 

6e^ 


(57) 


which is valid in the limit as e ^ 0 [82]. This calculation supports the hypothesis that the 
rate of caustic formation has a non-analytic behaviour as the Stokes number approaches 
zero. 
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